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Abstract 

(~ i I In this paper we present two atomistic models for the energy of a one-dimensional 

CL(. elastic crystal. We assume that the macroscopic displacement equals the micro- 

scopic one. The energy of the first model is given by a two-body interaction po- 
tential, and we assume that the atoms follow a continuous and piecewise smooth 
Cy ■ macroscopic (continuum) deformation. We calculate the first terms of the Taylor 

expansion (with respect to the parameter representing the interatomic distance) 
of the atomistic energy, and obtain that the coefficients of that Taylor expansion 
represent, respectively, an elastic energy, a sharp-interface energy, and a smooth- 
interface energy. The second atomistic model is a variant of the first one, and its 
Taylor expansion predicts, in addition, a new term that accounts for the repulsion 
force between two sharp interfaces. 
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1 Introduction 



j-y-j ' The aim of this paper is to derive the continuum expression of the elastic energy and 

f^*) , interfacial energy of a one-dimensional elastic crystal from an atomistic model. The 

motivation of that analysis came from our desire to justify the continuum model proposed 

in Ball & Mora-Corral [3], according to which the same material can exhibit smooth 

j^ ' and sharp interfaces. In the one-dimensional case, that model is briefly described as 

H I follows. The elastic solid is represented by the interval (a, b), for some a < b. An elastic 

deformation of the body (a, b) is represented by an increasing, absolutely continuous map 

u : (a, 6) — >• K such that the function u restricted to (a, b) \ S is in the Sobolev space 

W^'^, for some finite set S depending on u. Then, in [J] we postulated that equihbrium 

configurations are minimisers of the energy / defined by 

I{u):^ ( [W{u'{x))+e'^u"{xf]dx + KCa.T:A{Su'), (1) 



where Su' is the set of discontinuity points of u', and e, k > are two small parameters. 

In p|l we saw that, when W is assumed to have two wells, ifO<K^e^l then the 

global minimisers u of / (subject to appropriate Dirichlet boundary conditions) satisfy 

Card 5*11' = 1, and, in particular, they present sharp interfaces, whereas if < e <C k <C 1, 
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then the global minimisers satisfy Su' — and present smooth interfaces. This suggests 
that in the expression ([1]), one should replace k with a multiple of e, since only when k 
is comparable to e can both smooth and sharp interfaces appear simultaneously. 

In this paper, we try to derive the energy ([T]) (once k has been substituted by a 
multiple of e) from an atomistic one. The main idea, which is based on the procedure 
of Blanc, Le Bris & Lions [6], is the following. We assume that there is a continuum 
deformation w : [a, 6] ^ M such that, for every small £ > 0, the atomistic (discrete) 
deformation is the restriction of u to [a, b] n eZ, and the atomistic energy of that discrete 
deformation is given by 

^ -' i#iezni[a,fc] ^ ^ 

where W : K\ {0} ^ R is a two-body interaction potential with suitable decay at infinity 
(for example, a Lennard-Jones potential). Next, we do a Taylor expansion of E^{u) with 
respect to e, with the hope of recovering the expression ([1]). To be precise, we compute 
the first terms of that expansion: 

Ee{u) = Ea+Eie + E2e'^ +o{e'^), (3) 

and, in this way, we partially succeed in justifying the model ([1]), in the sense that Eq 
accounts for the elastic energy, Ei accounts for the sharp interface energy (plus boundary 
terms), and E2 accounts for the smooth interface energy (plus boundary and jump terms). 
The essential difference with the paper of Blanc, Le Bris & Lions Ji^ is the regularity 
assumption on u. While in [B] they assumed that u was smooth enough so that all the 
Taylor expansions made sense, here we assume that u is continuous and piecewise smooth, 
so that u' can have jumps. In this way, we obtain, within _Ei, a term that accounts for 
the sharp interfaces. Sections [2]-l6] are devoted precisely to the computation of ^ from 

In Sections [7] and [5] we change slightly the above atomistic model in order to get, in 
the continuum limit, an additional term that accounts for the repulsion between sharp 
interfaces. The motivation of that analysis was led by our desire to obtain the expres- 
sion of a repulsion term between interfaces from an atomistic energy. In the continuum 
context, it is not clear how to define an energy that represents a repulsion term between 
interfaces, let alone in our setting where we have two kind of interfaces: smooth and 
sharp. In general terms, that repulsion energy would be given by a non-local expression 
depending on the distance between interfaces, and be a decreasing function of that dis- 
tance which tends to infinity when the distance tends to zero. It is not clear, however, 
how to express that in a formula, since a 'smooth interface region' is only defined in 
vague terms as a region where the absolute value of the second derivative is very high. 
Among the three possible kinds of interaction between two interfaces (smooth-smooth, 
smooth-sharp and sharp-sharp), only the sharp-sharp interface interaction energy is easy 
to model, namely, as a function depending on the distance between the two interfaces, 
as describe above. Here we are using the fact that, in dimension 1, a sharp interface is 
represented by a single point; in higher dimensions, in contrast, the sharp-sharp interface 
interaction energy is not so easy to define. 
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In Section [7] we try to derive this repulsion term between sharp interfaces, and, 
again, we succeed only partially. The derivation follows a similar procedure to the one 
described above, but with an important difference that we describe now. Let e > 
be the interatomic distance. Then we assume that the two sharp interfaces are located 
at to G (a, 6) and io + "i£ G (a, 6), where jn > 2 is a natural number. So we are 
assuming that the interfaces are at a distance which is a fixed multiple of the interatomic 
distance. This assumption was motivated by the experiments of Baele, van Tendeloo & 
Amelinckx [2] , who observed a quasiperiodic microtwinning in the alloy Ni-Mn resulting 
from a Martensitic transformation, and whose images suggest that two consecutive sharp 
interfaces are separated by a distance which is 6 or 10 or 11 times the interatomic distance. 

The procedure is as follows. As above, we assume that there is a continuous deforma- 
tion u : [a,b] —>-R which is smooth in [a, to] and [to, 6], for some to € (a, b), and such that 
the atomistic deformation Ug follows u in [a, to] U [to + sm, b] (much like in the analysis 
of the Section 11]), but follows a discrete deformation y : {to + e, . . . , to + (m — l)e} -^ M 
in (to, to + me). The discrete deformation y will be given by an optimal profile problem. 
To be precise, we consider the energy ^ of this Ue and calculate its Taylor expansion 
([3]). We interpret the difference between this energy and the energy of the first model 
as a repulsion term between interfaces. Naturally, the coefficients Eq, Ei, E2 will depend 
not only on u but also on y. We will see that the coefficient Eq does not depend on 
y. The coefficient Ei does depend on y, and we choose y to minimise Ei. It turns out 
that, in many cases, the optimal y is the straight line. This might explain why, in the 
region between two sharp interfaces, the atoms are aligned in a straight line. A further 
minimisation process shows that the optimal m is 6, which qualitatively coincides with 
the experiments of Baele, van Tendeloo & Amelinckx [2] explained above. 

In this paragraph we compare our approaches to similar ones found in the literature. 
As mentioned above, the analysis of Section |4] follows closely that of Blanc, Le Bris & 
Lions j6j , with the important difference that our continuum deformation u is continuous 
and piecewise smooth (instead of being smooth). The closest to the analysis of Section [7] 
that we have found in the literature is the paper of Blanc & Le Bris [5] . They assume that 
the sharp interfaces are at a distance 7 which is larger than the atomistic scale but smaller 
than macroscopic, soO<e<^7<^l; this is different from our approach, as our 7 would 
be ms. Their atomistic energy equals the analogue of ([2]) plus a term accounting for 
the energy between the two sharp interfaces. That energy is given again by an optimal 
profile problem, but, in their case, the deformation between two sharp interfaces is a 
continuum one. Although they approach is very natural, it does not predict a repulsion 
term between interfaces; quite the opposite: the energy associated with two consecutive 
sharp interfaces is an increasing function of the distance between them. We believe that 
that conclusion is wrong, and this was one of our motivations to present our model of 
Section [T] 

Of course, there is a great number of atomic-to-continuum derivations for elastic and 
surface energy. We have pointed out those whose approach is similar to that adopted 
here. For radically different approaches, see, for example, Arndt & Griebel [1] (who use 
an upscaling method) and Braides & Cicalese 7| (who use F-convergence) , as well as the 
references therein. 

There are several disadvantages to be found in the atomic-to-continuum method used 
in this paper. Some of them are already present in Blanc, Le Bris & Lions [6]: the 
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macroscopic displacement is assumed to be equal to the microscopic one, the atomic 
deformation is assumed to follow a given smooth (or piecewise smooth, in our case) de- 
formation, the lattice is assumed to be periodic even near the boundary of the body, 
the Taylor expansion Q is not uniform in u, there is no guarantee that minimisers of 
the atomistic energy converge to minimisers of the continuum energy, and, finally, the 
limit continuum problem is not well-posed (because it is not coercive). In addition, in 
our approach we also have the disadvantage that, although we allow for more general 
deformations (piecewise smooth as opposed to smooth) , the discontinuity set of the de- 
formation gradient is prescribed. There are, nevertheless, some advantages to be found 
in our approach. First, it detects the right scaling between smooth and sharp interfaces: 
the term accounting for sharp interfaces appears in the coefficient e of the Taylor ex- 
pansion ([3]), while the term accounting for smooth interfaces appears in the coefficient 
e^; this corroborates the scaling deduced in Ball & Mora-Corral [4] by F-convergence 
methods. Second, it detects (or suggests) the 'general form' of a continuum energy func- 
tional accounting for elastic energy, sharp-interface energy and smooth-interface energy, 
thus partially justifying the model of 0]. Third, it is able to predict a repulsion term 
between sharp intrefaces and to describe roughly the atomistic configuration between 
two consecutive sharp interfaces, in qualitative agreement with the experiments of Baele, 
van Tendeloo & Amelinckx [2] . We do not consider the restriction to the one-dimensional 
setting to be a disadvantage, because we believe that many of the calculations in this 
paper can be carried out to the higher-dimensional case, as done in Blanc, Le Bris & 
Lions [5]. 

The outline of this paper is as follows. In Scction[2]we introduce the general notation 
of the paper, and state one of the main results of Blanc, Le Bris & Lions [6], which is 
the Taylor expansion ([3|) for smooth deformations, and constitutes the starting point of 
this paper. In Section [3] we study the map that transforms the atomistic potential into 
the elastic one. To be precise, a consequence of [6 (which we recall in Section [5]) is that 
the term Eq of ^ has the form J Wi{u'), for a function Wi depending on W. We show 
that blow-up rates at and decay rates at cxd of PF imply the corresponding properties 
on Wi, and we give sufficient conditions for the map W *—>■ Wi to be an isomorphism. 
Section m is the core of the paper and contains the calculation of the Taylor expansion ([3|) 
for a continuous and piecewise smooth deformation u having exactly one sharp interface. 
In Section [5] we compare the conclusion of the result of \Q^ recalled in Section [2] with 
our result of Section |4l In our analysis, there is a new term accounting for the sharp 
interfaces; we study the sign of that term and give sufficient conditions for it to be 
positive; this positivity physically means that we need energy to create a sharp interface. 
Section [6] shows that if the deformation has several sharp interfaces that are separated 
at a macroscopic distance, then the analysis does not essentially differ from the case 
of one interface (showed in Section |4]) , since the term accounting for the interaction 
between sharp interfaces is of order o(e^). Section [7] contains the other principal result 
of the paper: we present the model for a deformation with two sharp interfaces that 
are separated at a distance multiple of the atomic one, and calculate the corresponding 
Taylor expansion. Section [5] compares the conclusion of the results of Sections H] and [T] 
The new term that appears in Section [7] is a term that accounts for a repulsion force 
between the two sharp interfaces; in some particular cases, we solve the optimal problem 
and, thus, describe the atomistic configuration between two consecutive sharp interfaces. 
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2 Smooth deformations: no sharp interfaces 

In this section we introduce the general setting and notation of the paper, and state 
the resuh on the Taylor expansion of the atomistic energy in the case of a smooth 
deformation. 

It is known that in the atomic-to-continuum analysis, the Taylor expansion as s goes 
to zero depends on the particular sequence of e going to zero (see, e.g., Blanc, Le Bris 
& Lions [6] or Braides & Cicalese [7]). Usually, there is a natural choice of sequence 
of e going to zero. In this paper, we calculate all the possible limits according to the 
particular sequence of £ going to zero. In the next two lemmas, we introduce the language 
and main properties in order to deal with particular sequences of e going to zero. 

Lemma 1 For every a,b,c £ M. and £ > satisfying a < b, let fci, /c2, N £ Z be defined 
by the condition 

-a-c < fci < -a- c+ 1, -b-c-1 < k2 < -b-c, N ^ k2 - h + 1. (4) 

£ £ £ £ 

Then N = Card£Z n [a, b] and, as e ^ 0'^ , 
1 1 1 



+ 0(£), £(c+fci)-a = 0(£), 6-£(c+fc2) = 0(£). 



iV£ b- a 
Now let £ ^ 0+ be a sequence such that there exist 

01,02,61, 62, ci,C2 e M (5) 

satisfying 

e{c + fci) - o = ai£ + 02£^ + o(£^), b~ e{c + fc2) = &ie + &2e^ + o(£^), 

Ci£ + C2£^ +o(£^); 



11 1 . . -2 , ,., (6) 



then 



N e b-a 



< oi < 1, if ai = I then 02 < 0, if ai = then 02 > 0, 

< &i < 1, if bi = l then 62 < 0, if bi ^ then 62 > 0, 

|ci| < Yi v7' */ "^i = 71 V? ^^^^ "^2 < 



(6-o)2' ' ' {b-aY ' - {b-af 

*/ ''I ~ 71 72 then C2 > 



Proof. We have 



(b -af {b- ay 



b — a 
KN < + 1, 



hence 

-1 1/^11 



<- ^--. < 



{b — a + e){b — a) e \N e b — aj {b — a — e){b — a) 
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Ifci = 


= Tb^ *^™ 










1/11 


1 


-cie 


< 




b — a 


where 


as if ci = -(jfip- then 










1/11 

£2 i iV £ ' 


1 


-cie 


> 




b — a 



(6 — a — e) (6 — a) 



1 



2 ' 



(6-a + e)(5-a)2' 

Similarly, 

— oi £(c+ fci) — a — flie 1 — ai 

and this implies the corresponding properties for ai,a2. The proof for bi and 62 is 
analogous. ■ 

For the sake of simplicity, many of the theorems in this paper will be stated and 
proved when the deformation is defined in [—1, 1], and the lattice is Z. For this choice, 
we can say more about the parameters that appear in Lemma [H 

Lemma 2 Let a = —1, 6 = 1 and c — 0. For each e > define ki, k2,N ^ Z by (|3]). 
Then fc2 = — ^i cind N ~ 2fc2 + 1. Now let e — > 0+ be a sequence such that there exist 
^ satisfying ©; then 

ai 1 a? ai a2 1 

b,=a,e[0,l], b,^a,, ci = ---, c,^--- + - + -, ^^^ 

if ai = 1 then 02 < 0, if ai ~ then 02 > 0. 

Moreover, for every (O such that ([7]), there exists a sequence £ — > 0+ such that ^. 

Proof. The equalities ^2 = — fci and N — 2^2 + 1 follow at once from the definition. 

Now take ([5]) and a sequence £ -^ 0+ such that ^. The facts 61 = ai and 62 — 0,2 
follow from the equality ^2 = — fci. We have, successively, 

£^2 = 1 - ai£ - a2£^ + o(£^), iV£ = 2 + (1 - 2ai)£ - 2a2£^ + o(£^), 

and using the general formula 

(po+Pi£+P2e'+o(£'))^i ^ L-Ple+PLZ^^e'^+oie''), po e M\{0}, pi,p2 e M, 

Po Po Po 

(8) 
we obtain the equalities of ci and C2 of ([7]). The rest of relations of ([7]) follows from 
Lemma [T] 

In order to prove the last part of the theorem, we construct, for every ai G [0, 1] and 
a2 G M such that 

if ai = 1 then 02 < 0, if ai = then 02 > 0, 
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the sequence {enjnGN of positive numbers tending to zero defined by 

/ „^+a"n+a. '^^ K , a^) ^ (f , 0) , 

I Tr^^^T-, if (ai,a2) = (l,0) 

The corresponding sequence {fc"}„gN of fci satisfies fc" — —n for big n G N, and £„fc" + l = 
aie„ + a2e„ + o(e^). ■ 

With all these preliminaries, we are in a position to present the model of the atomistic 
energy of a deformation of an elastic crystal. The rest of the paper will be, of course, 
devoted to the analysis of that model. 

The one-dimensional elastic body is represented, in the reference configuration, by 
the closed interval [a, 6], for some real numbers a < b. Although it is customary to 
represent the reference configuration through an open set, in our case, we believe that 
the calculations are slightly simpler; in any case, it makes very little difference. The 
continuum deformation of the body is represented by a continuous increasing map u : 
[a, b] — > M. The continuity models that no fracture is allowed, and being increasing models 
the orientation-preserving character of the deformation and the non-interpenetration of 
matter. We assume that the body possesses a crystalline structure; in particular, we 
choose the lattice £ := c + Z for some c G R, we take £ > as the interatomic distance 
(which in the end will go to zero) and assume that the atoms of the body are located 
at the points of ei n [a,b]. Thus, the number of atoms of the body is Card£^ n [a,b], 
which, in the notation of Lemma [H coincides with N. We assume that the atomistic 
deformation u^ : e£n[a, 6] ^ R is the restriction to eiri[a, b] of the continuum deformation 
u. The atomistic energy of the discrete deformation u^ is given by a two-body interaction 
potential. This assumption is known to be very simplistic, but we believe that a good 
understanding of this model is needed prior to the analysis of more general and realistic 
ones. So let W : M\{0} ^ R be the two-body interaction potential, which is a continuous 
function with some decay properties at infinity; the precise assumption will be stated in 
Theorem [3] below. We assume that the atomistic displacement equals the macroscopic 
one, and, hence, we define the atomistic energy of the deformation u^ as 

j^ I \ 1 >r- Tx. /^ Me(£j) -u^(ei) \ 
^-("-^•^ 2Card£^n[a,6] \ ^ (, S j" ^'^ 

Note that, in terms of the notation ([3]), the energy can be equivalently written as 

^ J_ ''y^' / u(£(c+ fci -I- j)) - m(£(c+ fci + i)) 

The above paragraph has described, essentially, the particularisation to dimension 1 
of the model presented by Blanc, Le Bris & Lions [B]. The only modification is that, 
in our case, the body is represented by a closed interval (not an open one) and we are 
working with an arbitrary sequence of £ ^ 0+. Because of those (minor) modifications, 
the following result is not a particular case of Blanc, Le Bris & Lions pi Th. 3], but, 
since the proof follows exactly the same lines (and in fact, it is simpler because of the 
ID assumption), we omit it. 
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Theorem 3 Let W : R\ {0} ^ M. be a C°° function such that W{x) = W{-x) for all 
a; e M \ {0}, and satisfy that there exist C, R > and a > 3 such that 



\W''> {x)\ < C\x\-°'-\ for all xeR\{~R,R) and i £ {0,1,2, .. .}. 



(10) 



Let a,b,c € M and e > with a < b. Define £ := c + "E. Let u : [a,5] — > M 6e a C°° 
difjeomorphism. Let Ue he the restriction of u to eH. r\[a, b] Take ^ and a sequence 
e -^ 0+ such that ^. Define ©. Then 



where 



E^(ue) = E"+ eE^ + e^E^ + o{e^), 
i?":=- / J2w{u' {x)j)dx, 



(11) 

(12) 



E^ := ci{b - a)E" - 



1 



2(5 -a) 



J2 [iJ + 2ai - l)Wiu'{a)j) + (j + 26i - l)W{u'{b)j)] , 



j=i 



E' :=- 



1 r^ °° 

^i(63^ / E W"{u\x)j)u"{x)'f dx + C2{b - a)Eo 



(13) 



E 



fliCi 



Eiy-^^^^ 



i=i 
1 
b — a 



0-2 

b — a 
fe- 



ci 



J W{u'{a)j) 



- 'ij) Wiu'ib)j) 
a 2 / 



(14) 



b- 



— E 

1 °° 

— E 



J = l 



1 

"l2 


+ 


ai 
2 


- 


1 
L2 " 


2 


+ 


5? 
2 



J + 



j + 



2 

&i 
2 






M^'(«'(a)j>"(a) 
M^'(7.'(5)j>"(&). 



Theorem [3] requires some regularity of u and W , and decay conditions on W and its 
derivatives. As explained in the proof of Blanc, Le Bris & Lions [6, Th. 3], it is possible 
to prove pTjl with a different set of hypotheses on u and W; roughly speaking, less 
regularity of u requires stronger decay conditions on W . In fact, a proof of a version of 
Theorem [3] with weaker assumptions on u and W was made by Blanc & Le Bris [5, Th. 
2.1]. However, in this paper we are not interested in relaxing the regularity assumptions 
on u. Note also that, in Theorem [31 we have written explicitly the boundary term in 
P^ . which was not done in the higher-dimensional case of Blanc, Le Bris & Lions ^. 

3 The map that transforms the atomistic into the 
elastic potential 

In Nonlinear Elasticity Theory (see, e.g., Ball [3j), it is assumed the existence of a 
function W\ : M ^ K U {oo}, called the (elastic) stored-energy function of the material, 
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such that the elastic energy of a deformation u : (a, 5) ^ R of the body represented by 
the interval (a, b) is / Wi{u' {x))dx (of course, we have restricted the general theory to 
the one-dimensional case). This and Theorem [3] suggest the introduction of the operator 
that maps any function W : (0, oo) -^ M into the function Wi defined by Wi{t) := 
J27Li W^(i^)i for each t > for which the series converges. In this section we give a 
sufficient condition for that operator transforming the atomistic potential W into the 
(continuum) elastic one Wi to be an isomorphism. 

Let p,q eM.. Define Ap^q as the set of / £ C(0, oo) such that limsup(_,o+ ^''1/(^)1 < oo 
and limsupj^QQ *'' 1/(^)1 < oo. Clearly, Ap^q is a vector space. For each / £ Ap^q, define 



\f\\p,. 



max ■ 



sup tP\f{t)\, sup i«|/(i)| 



te(o,i) 



te[i,oo) 



It is immediate to see that \\ ■ \\p^q is a norm in Ap^q that equips it with the structure of 
Banach space. 

In this section (and also in Section [5]), we let C : (l,oo) -^ M. denote the restriction 
to (l,oo) of Riemann's zeta function, i.e., C(s) '■— J27LiJ~^ ^or each s > 1. We denote 
the norm of a linear operator between two Banach spaces simply as || • ||; the identity 
operator is denoted by I. Motivated by the introduction of this section, for each fc G N we 
consider the operator Tk that maps any function / : (0, oo) — > KU {oo} into the function 
Tkf defined by Tkf{t) := X^^i j'^fiJ^) for each t > for which that series converges (to 
a number or to oo). 



Lemma 4 Let fc G N and p,q > k + 1. Then Tk 
operator. If p > q then \\Tk — I\\ < ({q — fc) — 1. 



A. 



P:9 



A 



p,q 



linear bounded 



Proof. Call Sk := Tk - /. For each t > 1 and / e Ap.q we have i« |S'fe/(t)| < 
{({q — fc) — 1) ||/||p,g. For each < f < 1, define jt as the only integer satisfying jtt > 1 
and {jt - l)t < 1. CaU At := E -'="2' j"-'" + tP~'^ Y.T=j, j'"' and A := sup,g(o,i) At. Then 
i^ \Skf{t)\ < ^t||/||p,g- If we prove that ^ < 00, then we will have shown that Sk is a 
linear bounded operator with ||S'fc|| < max{(^((j — k) — \,A}. 

Now we prove that A < 00. li p> q then A < ({q — fc) — 1, whereas if p < g then 



sup At < sup 

te[i,i) tG[i,i) 



j=2 j=jt 



<29-P(C(p-fc)-l) 



and 



/■oo 

sup At < (ip -- k) - 1 + sup tP-'i / s''-'' ds < C(P - fc) - 1 



2q-p 



te(o,i) 



q-k-1 



Therefore, A < 00. m 

Let fc e N. For each i G {0, . . . , fc} let pi, qi G M, and define p := {po, . . . ,pk) and 
q :— {qo, ■ . . ,qk)- Let Ap^q-k be the Banach space of functions / S C'^(0, 00) such that 
/') e Ap^^q, for each i e {0, . . . , fc}, equipped with the norm ||/||p,9;fe := X^Lo ll/*''llpi,9i- 
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Proposition 5 Let k E N. For each i G {Q,...,k}, consider pi,qi > z + 1. Define 
p := (pq, . . . ,pk) and q := {qo, . . . , qk)- Then Tq : Ap,q-k -^ -^p,g;fe ^^ o, hounded linear 
operator. If Pi > qt > (^^{2) + i for all i £ {0, . . . , k}, then Tq : Ap^g-^k -^ "4p,9;fc ^■s 'J'^ 
isomorphism. 

Proof. Let / e Ap^q-k- By Lemma HI 

\\Tof\U,k = ^ WiToffWp,,,, = E WTJ'^Wp^.'i^ < E WTMf^Wp.,,. < max m\\\\f\u,k. 

i=0 4=0 i=0 

If Pi > q.i > C"^(2) + i for all i £ {0, . . . , fc}, then, again by Lemma gl HTq - /|1 < 
maxo<i<fe \\Ti — I\\ < maxo<i<fe CiQi ^ *) ~ 1 < Ij and, hence. To is an isomorphism. ■ 

Note that C^^(2) — 1.72865. As an example, motivated by the Lennard- Jones poten- 
tial (see (1^ below), for each k G N define LJk ■= A(^i2^...^i2+k),{6,...,6+k);k- Then, the 
Lennard- Jones potential (|^ belongs to LJk, and by Proposition [5l Tq : LJk -^ LJk is 
an isomorphism. 

Finally, we recall that Ventevogel [9] constructed an example of a continuous function 
(/) : (0, oo) ^ M U {oo} with exactly one relative minimum, and such that TqcJ) has several 
relative minima. In fact, his example can be easily adapted to construct, for eachp, q > 1, 
a smooth function W 6 Ap^q such that W has exactly one relative minimum, and TqW 
has several relative minima. 



4 Piecewise smooth deformations: sharp interfaces 

This section is devoted to the proof of the following result, which is the analogue of 
Theorem |3] for deformations that are continuous and piecewise smooth. 

Theorem 6 Let W : R\ {0} ^ R be a C°° function such that W{x) = W{-x) for all 
X CzM\ {0}, and satisfy that there exist C, R > and a > 3 such that (jlOp . Fix a = — 1 
and b = I. Define i :— Z. Let u : [—1, 1] ^ M &e continuous, increasing and satisfy that 
w|[-i.o] o.'^^d u|[o.i] are C°° difjeomorphisms. For each e > 0, let u^ be the restriction of 
u to e£ n [a,b] Take ^ and a sequence e ^ 0~^ such that ([5]). Define ^. Then (jlip . 
where 

«1 oo 



1 f^ 

•^ — 1 „■ 1 



'-^■j=i 



E' := L -1)e"-\ E(^' + 2«i - 1) {W{u'i-m + Wiu'im] 

- EO- - 1) [w{u'{o-)j) + w{u'{o+)j)] + ^ E ^ ("'(0^)^' + "'(0^)* 



j=2 i,J = l 
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+ J E [-7 + «i - «2 - «i + (J - y) jI [W{u'{-m + W{u' {l)j)] 

+ 1 E [(^ - «i + «?)^' + («i - J'' + 1^'] [w'{u'{i)j)u"{i) - w'iu\-i)jy'i-i)] 

Y^ij -l)[W{u'{0-)j) + W{u'{0+)j)] 
2 E (p - 1^'' + J^'') [W'{u\0-)j)u"{0-) W\u'{0+)j)u"{0+)] 



i=2 
2 4 



1 \ °° 



*j'^l 



- 5] K'(0+)j2 - i."(0-)*2] VF' {u'{0+)j + u'{0- 



»J=1 

Proof. In ([4]), given a, 6, c G R and e > with a < 6, we defined fci, fc2 and A^. In our 
case, we have a = —1, b — 1 and c = 0, so define fci, fc2, iV accordingly. Define a" = —1, 
6^ = and c^ = 0, and construct k^ , fcj", iV~ accordingly; finally, define a+ = 0, 5+ = 1 
and c+ = 0, and construct k^ , A:^, N^ accordingly. It is easy to see that k^ — ki, 
k^ ^ 0, k^ ^ and k^ = k2, hence N^ ^ N - k2 and 7V+ = TV + ki. By Lemma El 
fc2 = -fci, TV = 2/c2 + 1, TV" = A^+ = ^ and ^. Define 

flj^ := ai, flj" := a2, 6]^ := 0, 63^ := 0, a^ := 0, a^ := 0, b^ := 61, 



6+1=62, Ci:=2ci--, C]^ := c^ , C2 := -2ci + 2c2 + -, cj := C2 . 
Then, straightforward calculations using Lemma [2] and Formula ^ show that 



(15) 



£(c± + kf) - a± == af e + a^e' + o(£^), 6± - e(c± + fcf ) = 6f e + 6f e" + o(e^) 



1 1 _ 1 

7V±e ^ 6±-a± 



+ cf£ + c*£^ + o(e^) 



Define 






P±. ^ 1 " V^ „. (u{s{kf+j)) - u{e(kf + z)) 



2iV± 






By Theorem[2]and dH]), £;±(u) = (£;°)± + e(^i)+ + e2(£;2)± + o(e2), where 

/O CO „1 00 

^W^K(a;)j)da;, (i?")+ := / ^ W^K(a;)j) da;. 
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_. OO _j oo 

{E^ := (ai - l){E^r - ^J^ij + 2«i - l}W{u' {-l)j) - ^Y.^] - l)W{u'{Q-)j), 

-l OO -.00 

(i?i)+ := (ai - l)(i?°)+ - - 5] (.7 - l)W^K(0+)j) - - 5] (.7 + 2ai - l)W^(«'(l)j), 

"'-1 ,=1 

W{u'{-l)j) 

1 oi /l ai 



(i?^)- 



24 

00 

E 



1 3ai 2 , /^l «i\ ■ 



+E 

00 

+E 



2 + y + U^y'^ 



w^(w'(o-)j) 



.i_,l4),,(^_l).,4,3 



W-'(u'(-l)j>"(-l) 



T.[h- V + ^^'^ ) w^'K(o-)j>"(o-), 



i=i 



12" 4" 6" 



(i?2)+ := - - / Y. M^"("'(2^)J>"(^) Y dx + (1 - 2ai + a^ + a?) (ii;") 

VKK(0+)i) 



1 a\ ( \ a\ 



E 

OC 

E f-4?^ + -f \f' ) VK'(z.'(0+b>"(0^ 



1 3ai 2 , /^l ai\ • 

22 ^V2 2/-^ 



00 

E 



12" 4" 6 

1 fli a? \ /I ai \ 9 1 o 

4 2/-^ 6-^ 



12 2 2 '-^ 



W-'K(l)j)w"(l). 



We express 
1 



Ee{u,) 



Clearly, 



2iV 



— A;i k2—ki —ki k2—k\ 

E+ E +25: E 

i,j—0 i.j — — ki-\-l i— Oj — — fci + l 



-fcl 



VK 



u(e(fci + j)) -u{eiki +i)) 



(16) 



1 g ^ O^M^i±il)-^M^i±^ . ^^,(.), 



2iV 






2N 



Continuum limits of atomistic energies: smooth and sharp interfaces 



13 



and, by ®, and ©, 



iV + 1 1 1 fai 1\ 2 / 2n 



2N 2 + 4^+ T-s'^ +°(^ 



Therefore, 



^2.^1^ =2^^^ +' 

i,j=0 ^ ^ 



hE'Y^hE^r 



ii^i 



(17) 



■ £ 



"' ^^{E^ + liEr + liET 



+ oie'). 



Now we observe that 
2N ^ [ e 



^ + \Etiu^)^^±W«'^^-<'^ 



2N 



3 = i 



with 



2N " ^ ' 2^ ' 



^-iET + \iEy 



'i-l)(^'r + yr + li^'y 



+ o{s'). 



By Taylor expansion, for each j G {1, . . . , fc2}. 



W 



u{ej) - u{0) 



W{u'{0+)j) + W'iu'{0+)j) 



u{ej) - u(0) 



-u'iO+)j +Oie')j 



2 \, -4-02 



W{u'iO+)j) + e-W'iu'iO+)j)u"{0+)f + 0{e') {f-^- + j'-'^^) 



Therefore, 

N "^ 



It is easy to see that 






i=i 



fe2 



^ W{u'{0+)]) = ^ M/(«'(0+)j)+o(£), ^ W'{u\Q+)j)f = ^ W-'K(0+)j)j2+o(l) 
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Thus, 



l.±w('-!i^l^)^eltwiu'(0-)j) 



J = l 



i=i 



Ol 1 

y " 4 



oo _, oo 

J2 W{u' {0+)j) + -J2 W\u'{0+)j)u"{0+)f 



j=i 



j=i 



+ o(e2). 



In total, 



k2 —ki 



— y^ w 



2iV 



ij = -fci+l 



i(i?V + £ 



u{e{ki +j)) -u{e{ki + i)) 



1 1 ^ °° 



j=i 



(18) 



'i-l)i^T^\iEr^liE^rU-'i^\ 



(et + -AE'r + kE^r + -T + 1 E w{u'{o+)j) 



J=l 



lf;Ty'K(0+)j>"(0+)/ 



j=i 



+ o(£2). 



Take < i < -fci and -fci + 1 < j < fca - fci. Then 

u{e{ki+j)) -u{e{ki+i)) z.^+w, , .^ , ,.„-w , .^ 
= u (0^)(fci + j) + ?i (0 )(-fci - z) 

+ f K'(0+)(fci + jf - u"(0-)(fci + if] + 0(e2) [(fci + jf + (-fci - z)3 



Therefore, 

^, u(£(A:i+j)) -u(£(/ci+z)) 



") = M^ {u'{Q+){ki + j) + «'(0")(-fci - ^)) 



+ w {u'{o+){k, + j) + w'(o-)(-fci - 1)) ^ K(o+)(fci + jf - u"(o-)(fci + if] 

+ 0{e^) ((.? - z)-"^^ [(fci + jf + (-fci - if] + (j - z)-'^^ [(fci + jf + (-fci - iff) 
Now 

^E e' ((•?-*)""' [(^1 +^f + (-fci -*)'] + (.? -0""^ [(fci + J? + (-fci -*)']' 



iV- 



i=0 j = -fci+l 



TV 



^^' E E [(^' + ^r" (/ + ^^) + (.? + r "^ (/ + ^ Y] 
^ ^is' E E [(j + *)'""' + (■?■ + »)'""'] ^ ^2£' 

1=0 j = l 



l + ^(i4-"l+i5-"^ 



= o(£'), 
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for some constants Ci, C2 > depending on W, u, but not on i, j, e. It is easy to see that 



00 00 



— ki k^—k\ 

Y, Y, T4/(«'(0+)(fci+,7)+i.'(0-)(-fci-.))-^^iy(«'(0+)j+,i'(0-).)+o(e) 

i=Oj = -fci+l i=0 i=l 



Therefore, 



1 



— ki k2—ki 



nT.11 W{u'{Q+){k,+J)+u'{0-){-k,-^)) 

i=0 j = -ki+l 



1 /fli 1 I ^ 



Similarly, it is also easy to see that 

-fci k2 — ki 



i=0 j=l 



^ E E ^' ("'(0+)(^i + 3) + "'(0-)(-fci - ^)) K(0+)(fci + 3? u"{0-)ih + ^Y 

^ 00 00 

= ^ I E E ^' K(0^)J' + "'(0^)*) K'(0+)J' - u"{Q-)^'] + o{e^). 



i=0 3 = 1 
In total, 



^ 00 00 

I E E ^' K(o^)^' + "'(0^)*) K(o+)j' - ""(0^)^'] 



1=0 j=l 



+ o{e^). 



(19) 



Equalities HI]), (H?]), dUl) and ^^ conclude the proof. 



5 Sign of the jump term 

In this section we compare the conclusions of Theorems |3] and |6l The main difference 
in the assumptions is that in Theorem [3] only smooth deformations are allowed, while 
in Theorem [5] we allow deformations that are continuous and piecewise smooth. In 
the corresponding Taylor expansion (jlip , under the assumptions of Theorem [5] this is 
reflected in the appearance of a jump-derivative term in the coefficient of order e and 
higher. As explained in Section [l] this term models the sharp-interface energy. For 
physical reasons, we believe that this term should be positive, so that we need energy 



16 Mora- Corral 



to create a (sharp) interface. In this section, we analyse the sign of that term and show 
that, in many examples, it is indeed positive. 

Let W : M \ {0} — > M. and w : [— 1, 1] ^ M satisfy the same assumptions of TheoremlS] 
Let _E°, E^, E^ be the coefficients defined in Theorem[6l Let E°, E^ , E^ be, respectively, 
the coefficients ([T2|) . ([T3|) . (fTl]). once the substitutions 



fli 1 a? fli a2 1 

a = -l, 5=1, 5i=ai, 63 = 02, ^i = y - |, ^2 = y - y + y + ^ 

have been made. Of course, those substitutions are motivated by Lemma [2] For each 
i e {0, 1, 2}, call J' := E' - E\ Then jO = 0, 

J'^--JZ^3-l)W(u'{Q-)J)—J2{3-^)W(u'{Q+)3) + -Y,W{u'{Q+)J+u'{Q-)^), 

-^^ =^ E (i - y) (J' - 1) [w^K(0")j) + w{u'{Q^)3)] 

j=2 ^ ^ 

+ ^E (^^' - i^'' + \^ [M^'(«'(0-)j>"(0-) - M/'(u'(0+b>"(0+)] 

^ 00 

+ i E K(0+)^'' - ""(0-)*'] W^' ("'(0+)j + ^.'(0-)^) . 

The equality J" = expresses the fact that jumps in the derivative do not affect the elastic 
energy, which corroborates the model ([1]), and is a known result in the F-convergence 
approach of the problem (see, e.g., Braides & Cicalese [7]). Thus, J^ seems to represent 
(an approximation of a scaling of) the sharp-interface energy. 

In order to ascertain the sign of J^ we define J : (0, 00)^ ^ M as 

-. oc -.00 -.00 

J{a,h):=—J^{3-\)W{aj)--J^{3-\)W{h3)^-Y^W(h3+ai), a,b>0. (20) 

J=2 j=2 ij = l 

We believe that nothing can be said in general about the sign of J, so we restrict ourselves 
to the analysis of the case when the jump of the derivative is small, i.e., when a ~ 6. 

Proposition 7 Let W G C^((0, 00)) satisfy that there are C,R > and a > 2 such that 

W'\t) <Ct-°'-\ t>R, ie{0,l,2}. 
Define J : (0, oo)^ -^ R as ^, and ^ : (0, 00) -> R as 

Aia):=j^Y.i3-3')W"ia3), a>0. (21) 

J=2 
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Then there exists a neighbourhood U of {{a, a) G (0, oo)^ : A{a) > 0} such that J{a, 6) > 
for all (a, b) G lA. 



2 ji„ „^ _ I A{a) -A{a) 



Proof. Some easy but tedious calculations show that, for all a > 0, 

J{a, a) = 0, DJ{a, a) = (0, 0), D\j{a, a) = , ^^^^ ^^^^ 

Elementary calculus then shows that the function G : (0, oo)^ -^ M defined as 

Gia b) ■= / S^ ^^ a, 6 > and a 7^ 6 
^ ' ' ' \ A{a) if a > 

is continuous. Hence there exists a neighbourhood U of {(a, a) G (0, cxd)^ : A{a) > 0} 
such that G{a, 6) > and J{a, b) > for aU (a, b) eU. m 

We finish this section with an example of an interesting potential for which the func- 
tion (pi]) can be computed. For a > 0, let W^ ■ (0, oo) ^ R be the Lennard- Jones 
potential defined by 

W^.(i):-(f)''-(9'' ^>0- (22) 

Some calculations show that if 

- > ., ' — TT^ ^ 0.603431 



a V 7[C(5)-C(7)] ; 

then Aa-{a) > 0, where Aa- is defined by (HH, but having replaced W with Wcr- Recall 
that the minimum of the Lennard-Jones potential is at 2^/^a ~ 1.12246cr, and, thus, in 
particular, Aa- is positive at that value. 

6 Several sharp interfaces, well separated from each 
other 

When the deformation u presents several sharp interfaces (i.e., u' is discontinuous at 
finitely many points), and the interfaces are well separated from each other (i.e., the 
points of discontinuity of u' do not depend on e), then we have an exact analogue of 
Theorem |6l and, in particular, in the Taylor expansion of the energy, there is no term 
accounting for the interaction between sharp interfaces. 

Proposition 8 Let W : R\ {0} ^ R be a C°^ function such that W{x) = W{~x) for all 
a; G M \ {0}, and satisfy that there exist C, R > and a > 3 such that (jlOp . Let n G N 
and to < ■ ■ ■ < ira+i. Define a := tg and b := i„+i. Let c G M and define i :— c + TL. Let 
u : [a, 6] — > R be continuous, increasing and such that u\u j. i is a C°° diffeomorphism 
for each p G {0, . . . , n}. Take ([5]) and a sequence e ^ 0^ such that ([6]) and, for each 
p€ {0,...,n}, 

e(c -I- fcP) -tp = a^e + a^e^ + o{e^), ip+i - s{c + fcf ) = b^e + b^e^ + o{e^), 
1 1 1 



NP e tp+i - tp 



+ cf£ + c^e2+o(£2), 
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for some af,a2,6^,62jC?jC^ ^ ^^ '^'^'^ where fc^,fcf,iVP are defined according to (g]), by 
replacing a with tp and b with ip+i. Let u^ be the restriction of u to ein [a, b], and define 
^. Then (dl]), where (O, 

p=0 ''*-P 3 = 1 

+ F2 («'(a), u\t^), u'itt), ..., u'{b), u"{a), u"{t^), u"{tt), ..., u"{b)) , 

for some dp G R depending on c\, c^, tp+i ~ tp, b ~ a (for p G {0, . . . ,n}), some 
Cp e R depending on ci, cf, C2, cf, ip+i — ip, fo — a (for p e {0, . . . ,n}j, and some 
Fi G C°°((0, cx))2"+2) anrf i^2 e C°°((0, cx))2"+2 x R2n+2) depending on W, a{, a\, b\, b^, 
c?, cf ; Wi - *p r/o^'P e {0, . . . , n};. 

Proof. The proof is very similar to that of Theorem [H] and will only be sketched. 
For each p G {0, . . . , n} and £ > define 

*#i6^ni[tp,tp+i] 

From Theorem [3] we know that Eliu) = {E'^f + e(£;i)P + e^[E?-y + o(e2), where the 
expression of {E!^f, (E^Y, (E^Y is given by ^, ^, dH]), respectively, but replacing 
a, b, ai, a2, 6i, 62, ci, C2 with ip, tp+i, a\, a^, b\, b^, c^, cf, respectively. We express 

n n— 1 

2/j-e^n^[a,fc] p-o 4^je^n^[ip,ip+i] p-o ze£n^[tp,tp+i] je^n^[ip+i,ip+2]\{i} 

n— 2 n 

+ 2E E E E 

p^o g-p+2 iemi[ip,ip+i] jG£ni[tg,ig+i] 

-E E E - E E • 

p=o ig<?ni[tp,tp+i] iG^n{ti,...,t„}\{i} »e£n{ti,. ..,*„} jet.n\[a.b\\{i} 
Now 

= ±'-^f^{Er + e±{'-^^{E^Y^ap{Er 

p=0 p=0 ^ 

+ -' E (H^ (^')' + "^ (^')' + '^^ (^°)') ' 

p=0 ^ ^ 
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where, for each p G {0, . . . , n}, the number ap G M depends on ci, c^, ip+i — tp, b — a, 
and the number /?p S R depends on ci, c^, C2, Cj, tp+i — tp, b — a. 
Arguing as in the proof of Theorem [6l it is easy to see that 



1 

m 



n — 1 n 

2E E E -E E E 

P=o ie£ni[tp,tp+i] jefni[tp+i,tp+2]\{*} p=o ie^ni[tp,tp+i] ie«n{ti,...,t„}\{i} 



' E E 

iGfn{ti,. ..,*„} je£ni[a,6]\{i}_ 
= B,{u'it^),u'{t+),...,u'{t-),u'{t+))s 

+ B2 {u'it^),u'{tt), . . . , u'(t+), w"(ir), ii"(i^), . . . , «"(4)) e' + o(e'), 

for some Bi e C°°((0, cx))^") and S2 £ C°°((0, oo)^" x M^") depending on W, a{, a^, &f , 
62, C]*, C2 for each p £ {0, . . . , n}, and ip+i — tp for each p E {1, . . . , n}. 
Finally, 

^'f i: E E H^( "'^^) ;"'") ), o(.--). (23) 

p=o ?=p+2 ie<?ni[tp,tp+i] ieini[t,,t,+i] 

and, of course, e"~-^ = o(£^)- ■ 

The reason why in Proposition [8] only local terms appear in the expansion of the 
energy but not an interaction term between interfaces is that the sharp interfaces are 
separated at a macroscopic distance. We will see in the next section that if the sharp 
interfaces are separated at a microscopic distance then a new term will appear in the 
Taylor expansion accounting for that interaction energy. 



7 Sharp interfaces separated at an atomic scale: re- 
pulsion term 

From Proposition [51 and especially from ([^5)) . we conclude that, in the model studied 
in Section [6] (which is virtually the same as that of Section |4]), sharp interfaces do not 
interact with each other. This is due to the decay conditions on W and to the fact that 
the sharp interfaces are separated from each other at a macroscopic distance. In this 
section, we will see that, if the sharp interfaces are separated at a distance comparable to 
the atomic one, then a small variant in the model predicts an interaction energy between 
two consecutive sharp interfaces. 

In this paragraph, we briefly explain the atomistic model of this section and how it 
differs from the model of Section [6] The assumptions on the potential W and on the 
atomistic energy Q are the same. As for the continuum deformation, we assume that u 
has two sharp interfaces, separated at a distance which is a multiple of the interatomic 
distance e. The main difference is that the atomistic deformation u^ does not follow u 
in the region between the two sharp interfaces. 
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Now we explain the model in more detail. We assume that the macroscopic deforma- 
tion u presents exactly two sharp interfaces, which are at a distance of a fixed multiple 
of the atomistic scale; to be precise, at a distance me, for some integer m > 2. It is not 
clear how to define the deformation in the region enclosed by the two interfaces, nor the 
energy associated to it. Here we assume that, to the left of the leftmost sharp interface 
and to the right of the rightmost sharp interface, the atomistic deformation u^ follows 
the continuum deformation u, which is everywhere continuous and of class C°° outside 0, 
whereas in the region between the two sharp interfaces, the atomistic deformation follows 
a scaling of a given deformation y, which in the end will solve an optimal profile problem. 
Thus, we assume that there exist an increasing homeomorphism u : [—1, 1] — > M such 
that u|r_i 0] and u|ro_i] are C°° diffeomorphisms, and an increasing function y : [0, 1] ^ M 
such that the atomistic deformation u^ : [— 1, 1] n eZ ^ R is defined by 

Ue|([-l,0]U[m6,l])n£Z = u|([-l,0]U[me,l])neZ, 

u{me) - u{0) , X u{0)y{l) - u{me)y{0) 

y(l)-y(0) me 2/(1) -y(0) 

Take such a y. Then, for all a > and 6 G M, the function ay + b gives rise to the same 
Ue- Therefore, we can assume, without loss of generality, that 2/(0) — and y{l) = 1. 
Moreover, y need not be defined in the whole [0, 1] but only on {0, —,..., 1}. Thus, 

j 
u^{je) — [u{me) — u{0)]y{ — ) + u(0), ^ < j < rn — 1. 

m 

The atomistic energy associated with the deformation Ue is still ([9]). Since the model 
differs from that of Section |4] only in what happens in (0, em), it makes sense to compute 
the difference of the energies between the two models. This is done in the next result. 

Theorem 9 Let VK : M \ {0} ^ M 6e a C°° function such that W{x) = W{-x) for 
all X E R\ {0}, and satisfy that there exist C,R > and a > 3 such that (jlOp . Let 
u : [—1,1] —>■ R be continuous, increasing and satisfy that u|[__i.o] and w|[o,i] o,i"e C°° 
diffeomorphisms. Let m > 2 be a natural number. Let y : { — , . . . , "'~^ } -^ (0, 1) be a 
strictly increasing function. For each e > 0, define u^ : [— 1, 1] D eZ ^ R by 

"£|([-l,0]u[me,l])neZ — ""I ([-l,0]u[T?xe,l])n£Z) 

j 

Ue{je) = \u{me) — u{0)]y{ — ) + u(0), 1 < j < m — 1, 

m 

and let 

^^^"^^ ■^ 2Card[-l,l]neZ ^ ^ [ e 

^^(") •=2Card[-l,l]neZ ^ ^ [ e 

Take ([5]) and a sequence e ^ 0+ such that ([6]). Then Es{ue)—E^{u) = eKi+e^ K2+o(s'^) , 
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where 



K 



^ oo m — 1 



i-0 j-1 

^ 771 — 2 777 — 1 



W u'{Q+)my{ — ) + u'{Q-)i - W{u'{Q+)j + u'(0")i) 



4=1 j=i+l 

771—1 CXD 



^ ^ w^U'(o+My(^)-2/(-)] 

Z ^ Z / \ m 777 



' m' 






7—1 j—m 



j ~ my{—) 
m 



-. oc 

-o("*-i)E^("'(o^)^')' 



J=l 



(24) 



OO 771 — 1 r 



'<''-\Y.Y. 



i=0 j=l 



W' (u'{Q+)my{ — ) + u'{0-)i 



u"{0+)m^y{^)~u"{0-)t^ 



777. — 2 771-1 



i—1 j^i-{-l 

777 — 1 OO 



m m 



(0^)E E W'iu\0+)m[y{^)-y{-)])m^[y{^)-y{-)] 



7TJ m 



J -my[ 



f - m^y{—) 
m 



i«"(o+) E E ^' h'(o 

i— 1 j^rn ^ 

iz."(0+)(m - !)£ W^'(u'(0+)j)j(m + j) + f ai - i) K, 



j=i 



Proof. As in Lemma [2l for each e > let /c2 be the maximum integer less than or equal 
to 1/e, and define N := Card[-1, 1] n eZ. Then A^ == 2fc2 + 1 and 



rn — 1 rn — 2 rn — 1 rn — 1 ^2 

E E + E E+EE 



1 

N 



i— — k2 :/ — 1 *— 1 :y— ■i+l "i— 1 j—fn 



^ u,{ej) - u,{ei) \ _ ^ f u{ej) - u{ei) 



(25) 



Take -fcz < i < and 1 < j < m - 1. Then 



^e(ej) - We(e«) u{em) - u(0) . j\ , u(O) - u{ei) 

= y — H 

£ erne 



u'{Q+)my{ — )~u'{0-)i+ „ 
m 2 



u"(0+)m22;(^)-u"(0-)i^ 
m 



+ 0{e^){-if 



and 



"^'^'^ "^"'^ = «'(0+)j - ^'(0-)* + '-[u"{Q+)f - u"{0-)^^] + 0(£2)(-z)^ 
e 2 
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Therefore, 

^ , Me(ej) -u,{ei) 



m 



+ 7;W ( u'{0+)my{ — ) - u'(0")i 



u"{0+)m^y{^)~u"{0-)f 



and 



W 



+ 0{e^) {{"if-"' + {-if-""') + Oie^)i-if-^'' 

u{£J) — u{£iy 



--W{u'{0+)j-u'{Q-)i) 



+ ^W'{u'{Q+)j - u'{0-)i)[u"{0+)f - u"{0-)i 



6 — 02 



Hence, using ^ and Lemma [H 

m-l r 



^ u,{ej) - u,{ei) \ _ ^ ( u[ej) - u{ei) 



i = -k-2 j=l 

oo m—1 



i=0 j = l 



W u'{0+)my{ — ) + u'{0-)i - W{u'{0+)j + u'{0-)i) 



2 4 ' 



OO m—1 r 

EE 

i=0 j=l 



W u'{0+)my{ — ) + u'{0-)i - W{u'{0+)j + u'(0")i) 



4=0 J = l 



W {u\Q+)my{ — ) + u'{0-)i 
m 



u"{0+)m^yi^)^u"{0-)f' 



m 



-W'{u'iO+)j + u'{0~)i)[u"{0+)f - u"{0-)i^]] + o(e2). 



(26) 



Now take 1 < i < m — 2 and i + l<j<m— 1. Then 

Ue{£j) ~ u^{ei) u{me)-u{Q) j i 

= y — - y[—)\ 

e e m m 



u'{0+)m[y{^) - y{-)] + u"{Q+)m'[y{^) - y{-)] + 0{e') 
m m Z mm 



and 



"^^^'^ "^''^ = u'{0+Kj - z) + '-u"{0+){f ^^) + 0(£2). 



Therefore, 

'"£(£j) - We(e«) 



W 



= W[u'{Q+)m[y{^)-y{-)] 
m m 



+ '-W' (u'iO+)m[yi^) - y(-)] ) u"iO+)m'[yi^) - y(-)] + 0(8^) 
2 \ mm/ mm 
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and 



W ( !^M_^iM') ^ Wiu'{0+){j - z)) + ^W'iu'{0+)ij - i))u"iO+){f - i") + Oie"). 



Hence 

_. m — 2 rn — 1 r 



N 



^ u,{ej) -u,{ei) \ _ ^ f u{ej) - u{ei) 



i—1 j—i-{-l 

■m — 2 m—1 r 

2—1 j—i-\-l 



W { u\0+)m[y{^) - y{-)] ] - W(u'{0+)(j - z)) 



m m 



^ N m — 2 m—1 r 

«1 _ i ^2 ^ j2 



2 4 



i—1 j—i-\-l 
m — 2 m—1 



W { u'{0+)m[y{^) - y{-)] ] - W{u'(0+)(j - z)) 
mm 



i—1 j—i-\-l 



W (u\0+)m[y{^) yi^)]j m^^) - v{^)] 
-W^'(u'(0+)(j-0)(j'-^')]+o(£^)- 



(27) 



Finally, take \ < i < m — 1 and m < j < k2. Then 



Uejsj) - Uejei) _ u{em) - m(0) , z , u{ej) - u{0) 
e s ' m e 



u'(0+) 



i - my{—) 
m 



+ f ""(0+) 



f -m^y{ — ) 
m 



2\,-3 



+ 0{e')j 



and 



£ 1 



Hence 



w^('^^iM^^^£M')^y^(,'(o+^ 



j - my{ — ) 



+ lW'[u\Q+) 



j -my{ — ) 
m 



u"(0^ 



f - m^y{—) 
m 



+ 0(e2)(/-"i +f'''') + 0{e^)f-^' 



and 
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Therefore, 



N 



m—1 k2 

EE 



^ ue{ej) -u^{ei) \ _ ^ fu{ej) - u{ei) 



2—1 j—in 

771 — 1 OC 



i—1 j—m 



W m'(0+) 



j - my( — ) 



W{u'{0+)U-i)) 



^ _ 1 ] p2 

2 4 ' 



EE 



i—1 j—m 

m—1 oo 



W u'{0+) 



j -my{ — ) 



no^) E E 

i—1 j—m 



W u'{0+) 



j — ■my{ — ) 



Wiu'iO+)U-i)) 



f ~m^y{ — ) 



(28) 



-W'{u'{Q+){j-z)){f-i^)\+o{e^). 
Equations ([251), ^^, ^^ and ^^ conclude the proof. ■ 



8 The optimal profile problem 

This section analyses the terms Ki and K2 of Theorem [3] As we explained in Section [71 
it is not clear how to define the atomistic deformation in the region between two sharp 
interfaces. In Theorem [9l we assumed that, in that region, the atomistic deformation 
Ue followed a scaling of a given discrete deformation y : {^, . . . , ^^^} — > (0, 1), but 
the actual values of y were left unspecified. We believe that y should be such that 
the difference of energy Eg{u^) — E^{u) is minimum. As this is a difficult problem, we 
approximate Ei;{ue) — E^{u) by its Taylor expansion, and, thanks to Theorem [SI we 
choose y to be a minimiser of Ki. For every m, this is a finite-dimensional minimisation 
problem, and it will turn out that, in many cases, the optimal choice of y is the identity 
map, which gives as the optimal value for Ki. 

As explained in Section [U the repulsion term between two sharp interfaces should 
be decreasing with respect to the distance between them, and tend to infinity as the 
interface goes to zero. Since we have here a discrete variable m that runs over {2,3,...}, 
the latter property makes no sense, but we still can expect that the repulsion term is 
decreasing with respect to m, at least for small values of to. In the context of Theorem 
[9l the repulsion term is defined as E^{u^) — Eg(u), but again we approximate it by its 
Taylor expansion eKi + £^K2. As explained in the previous paragraph, the term Ki 
provides us with information about the optimal shape of y (which in many cases turns 
out to be the identity), but does not give much information about the repulsion energy 
(because in many cases it is zero) or about m (because in many cases the optimal value 
of Ki is 0, regardless of m). Thus, we will study the term K2 and it will turn out that, 
in a particular but important case, the optimal m is 6. 

The rest of this section is devoted to making those ideas precise and giving some 
examples in which the terms Ki and K2 can be calculated. 

Given a function W : (0, 00) -^ R and a natural number to, > 2, we define 



i^m '-^ I (2^1 7 • ■ • I ^m-l) 



: < a;i < 



*\ Xf] 



<1} 



(29) 
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and Fm : (0, oo)^ x Km ^ M as 

oo rn— 1 

Fm(a, 6; xi, . . . , Xm-i) := X! X! [^ (fomXj + ai) - W{bj + ai)] 

(30) 

rn — 2 rn — 1 tti — 1 oo oc 

i— 1 j—i+1 i—1 j—m J — 1 

for each a,b > and (xi, . . . , Xm^i) G Wr?x for which all the series of ([50)) converge. Of 
course, the reason of this definition is that, according to Theorem [HI and specifically ([M]) . 
and following the notation there, 

K, = iF,„ ('„'(0-),z.'(0+);2/(l),...,2/(^^i^)) . 
2 \ m ra J 



p ?n — 1 



In the next lemma we study the minimisers of F^- As we will see, the point g™ e 
defined by 

g™ := —,■■•, (31) 

\m m J 

will play an important role. 

Lemma 10 Suppose that W G C^((0, oo)) satisfies 

lim W{t) = oo, limsup t" max{|W^(t)|, |VK'(t)|, |W^"(t)|} < oo, 

/or some a > 1. Let m > 2 be a natural number. Define (j29p . (j3ip and i^m : (0, oo)^ x 
Urn ^ K as (I30p . Let a, 6 > 0. T/ien f/iere exists a niinimiser of Fm{a,b; ■) in Um- 
Moreover, Fm (a, b; q,„) = and DFm (a, a; (jm) = 0. Finally, 

{oo 
2a2TO2 Vl^"(aj) i/ fc,£G {l,...,m-l} wii/i k = £ 

-a'^m^W" {a\k - e\) if fc, £ G {1, . . . , m - 1} with k ^ £. 

Proof. The assumptions imply that F,„ is of class C^, and, for each a, 6 > we have 
F„i{a, 6; x) — ^ oo as a; — > dhim with x G Z-/„i- This implies the existence of minimisers of 
Fm{a, 6; •) in Um- The rest of the lemma follows from a direct calculation. ■ 

The next step would be to compute the minimisers of Fm{a,b;-) and to ascertain 
whether or not the point qm is a (local or global) minimiser of F„i{a,a] •). The answer 
to these questions depends on W, a and 6; all we can say is that q„i is not, in general, a 
critical point of F„i{a, b; •) for a,b > with a ^ b. Since we do not think that there is a 
general answer to these questions, from now on we will concentrate on a specific example. 

Let cr > and let W^ : (0, oo) ^ M be the Lennard- Jones potential defined in (1^^ . 
From now on, for each natural tti > 2, let the function F^, defined in ([5(7| refer to the 
potential W„. Then, a direct computation shows that, for all a > 0, 

oo g _i4 g 
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1 2 
det D^i^3 ( a.a;-,- 



4ai2a-28 



,16 



4677752 



Scr^TT^ - 2079a^ 



9 (26c 



7a'' 



Based on those formulas and in Sylvester's Criterion, we find that 



D''F2 ( a, a; - j > if and only if ^ < (^^ 



1/6 



D'^Fs ( a, a; -, - I > if and only if — < 
3 3/ a 



Stt' 



467775 



78 



1/6 



225 



-21 



Here we are using the following notation: if yl is a symmetric matrix, by A > we 
mean that A is positive definite. Thus, we have a necessary and a sufficient condition 
for Qm to be a local minimiser of Fm{a,a; ■), for each m G {2,3}. In fact, numerical 
experiments with the software Mathematica [8] suggest that D^Fm{a,a;qm) > if and 
only if det D^Fm{a, a; qm) > 0, if and only if a < amiy, where the numerical value of «„ 
is shown in Table [T] 



1.2435 



1.2430 



1.2425 



1.2420 



m 


flm 


2 


1.24362 


3 


1.24280 


4 


1.24226 


5 


1.24192 


6 


1.24169 


7 


1.24153 


8 


1.24142 


9 


1.24133 


10 


1.24127 


11 


1.24122 


13 


1.24115 


15 


1.24111 


17 


1.24107 


19 


1.24105 



1.2415 



1.2410 



I I T T T T T 



10 



Table 1: Numerical values of i 



Recall that there are two natural values of a, namely, the minimiser of W^, which 
is 2^/^(7 ~ 1.12246(7, and the minimiser of the elastic energy (see Theorem [3]), i.e., of 

1382 \l/6. 



E^i W^^(j*): which is (j 



TTCT ~ 1.1193(7. Numerical experiments with 



Mathematica [8] suggest that qm is in fact a global minimiser of Fm(a(7, acr; ■), where 
5:=(glf§5)'/^,forme{2,3,4,5}. 

We sum up the findings up to now. In the model described in Section [71 the interac- 
tion energy term between the two sharp interfaces located at and me is, by definition, 
E^{u^) — Ei^{u), which depends on u, y, e, m (and, of course, W, but this is fixed before- 
hand). Theorem[5]enables us to approximate E^{u^) — E^{u) by eKi, so an approximation 
of a scaling of the interaction term is given by Xi, which depends on u'(O^), u'(O^), y, 
m. Now we minimise in y, and hence the interaction term depends on u'(0"), u'(O^), 
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m. According to the analysis above, we believe that in many interesting cases, qm is a 
minimiser of F„i{a,a;-). Assume that this is the case. Then, by Lemma [TOl the first 
order term Ki of the expansion of Theorem [5] is zero, regardless of to; hence, in order to 
have more information we study the second order term K2 ■ As before, very little can be 
said about that term unless we assume a specific form for the potential W. Motivated 
by the previous analysis, we define the function G : (0, 00)^ x {2, 3, . . .} -^ R as 



G{a, a, to) := ^ {-f + 2f - j) W^aj) - (to - 1) 



E 



{j~l){2j-m)W;^{aj), 



for each a, cr > and to G {2,3, . . .}. It is easily checked that the expression K2 of 
Theorem[9l when the quantities u'(0~), u'(0+), u"(0~), u"{0'^) have been replaced with 
a, a, p, p, respectively, and the functions W, y have been replaced with W^, id, becomes 
precisely pG{a, cr, to)/4. For each cr > 0, we choose Ua- '■— ( g^|Pg ) 
the values of G{aa,cr,m). It turns out that the quantity <TG{a„ 
on cr, and their numerical values for several m are displayed in Table [21 together with a 
graph of those values. Again, the numerical values, as well as some symbolic calculations 
involving the Riemann zeta function, were obtained using Mathematica ^. 



TTcr, and we compute 
, a, to) does not depend 



-0.045 h 



m 


aG{aa,<T, m) 


2 


-0.0570514 


3 


-0.0657517 


4 


-0.0470596 


5 


-0.0453827 


6 


-0.0452401 


7 


-0.0452798 


8 


-0.0453306 


9 


-0.0453703 


10 


-0.0453990 


11 


-0.0454194 


12 


-0.0454342 


13 


-0.0454451 


14 


-0.0454533 


15 


-0.0454594 


20 


-0.0454752 


25 


-0.0454809 


30 


-0.0454834 


40 


-0.0454854 


50 


-0.0454861 



-0.055 



-0.04525 
-0.04530 
-0.04535 
-0.04540 
-0.04545 
-0.04550 



10 



'TTTTTTTTTTTTTTTT 



20 



25 



30 



35 



40 



Table 2: Numerical values of (TG{a^, a, to.) 



Our interpretation is the following. For the Lcnnard- Jones potential, when u'(0 ) ~ 
u'(O^) ~ a/a ~ 1.1193, then the optimal discrete configuration between two sharp 
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interfaces separated at a distance me is close to qm- If, in addition u"{0~) ~ u"{0~^) > 
then the optimal m is 6, whereas if u"{0~) ~ u"(0"'") < then the models predicts that 
no interfaces at all (i.e., m = 0) is energetically better. So this model predicts that the 
optimal length of the space between interfaces is 6 times the atomistic distance, which 
coincides with the experiments of Baele, van Tendeloo and Amclinckx [2]. 
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